Scale dependence of cosmological backreaction 
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Due to the non-commutation of spatial averaging and temporal evolution, inhomogeneities and 
anisotropies (cosmic structures) influence the evolution of the averaged Universe via the cosmological 
backreaction mechanism. We study the backreaction effect as a function of averaging scale in a 
perturbative approach up to higher orders. We calculate the hierarchy of the critical scales, at 
which 10% effects show up from averaging at different orders. The dominant contribution comes 
from the averaged spatial curvature, observable up to scales of ~ 200 Mpc. The cosmic variance 
of the local Hubble rate is 10% (5%) for spherical regions of radius 40 (60) Mpc. We compare our 
result to the one from Newtonian cosmology and Hubble Space Telescope Key Project data. 
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Various cosmological observations, interpreted in the 
framework of spatially flat, homogeneous and isotropic 
cosmogonies, have now confirmed the accelerated expan- 
sion of the Universe. The most direct evidence comes 
from the study of supernova (SN) of type la Many 
attempts have been proposed to understand this mystery, 
e.g., dark energy in the form of a cosmological constant, 
quintessence field or modification of gravity. However, 
these suggestions always rely on the homogeneity and 
isotropy of the cosmic medium, which are rather rough 
approximations. 

The Universe hosts enormous structures. In our neigh- 
borhood, there seem to exist two voids, both 35 to 70 Mpc 
across, associated with the so-called velocity anomaly 
a large filament known as the Sloan great wall about 
400 Mpc long [3j and the Shapely supercluster with a 
core diameter of 40 Mpc at a distance of ~ 200 Mpc 
from us Q. Furthermore, based on the Hubble Space 
Telescope (HST) Key Project data @, evidence for a 
significant anisotropy in the local Hubble expansion at 
distances of ~ 100 Mpc was found [6|, and an anisotropy 
of SN la Hubble diagrams extending to larger distances 
has been reported recently Therefore, spatial ho- 
mogeneity and isotropy seem to be valid only on scales 
larger than ~ 100 Mpc @, and effects of local inhomo- 
geneities are worthy of investigation. More specifically, 
observables from within a few 100 Mpc must be revisited 
critically. The most fundamental of those are cosmic dis- 
tances and the Hubble constant Hq. 

In this paper, we study the averaging of the inhomo- 
geneous and anisotropic Universe over a local domain in 
space-time. We stick to the idea of cosmological inflation, 
assuming that the Universe approaches homogeneity and 
isotropy at scales as large as the Hubble distance. 

Many cosmological observables are averaged quanti- 
ties. For instance, the matter power spectrum is a Fourier 
transform and thus a volume average weighted by a factor 
e ik-x Another very important example is the idealized 
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measurement of Hq . One picks N standard candles in 
a local volume V (e.g., SN la in the Milky Way's neigh- 
borhood out to ~ 100 Mpc), measures their luminosity 
distances di and recession velocities = czi [zi being 
the redshift of each candle) and performs the average 
Ho = f 1 - I n the limit of a very big sample, it 

turns into a volume average Hq — y J gdV. 

Cosmological observations are made on the past light- 
cone, so one should average over a light-cone volume. 
However, for objects at z < 1, spatial averaging on a 
constant-time-hypersurface is a good approximation, as 
the Universe does not change significantly on the tempo- 
ral scale involved. 

Due to the nonlinearity of the Einstein equations, spa- 
tial averaging and temporal evolution do not commute. 
Hence, inhomogeneities and anisotropies affect the evo- 
lution of the averaged Universe via the so-called "back- 
reaction mechanism" [H El GJ, El EI M El- Below, 
we utilize Buchert's averaging method [12| to estimate 
the order of magnitude of backreaction effects and study 
the signatures of averaging from the local measurement 
of H . 

Buchert's setup is well adapted to the situation of a 
real observer, if we are allowed to neglect the difference 
between baryons and cold dark matter (CDM). On scales 
> 10 Mpc, baryon pressure is insignificant, and a real ob- 
server comoves with matter, uses her own clock and re- 
gards space to be time-orthogonal. These conditions de- 
fine a comoving synchronous coordinate system. There 
are no primordial vector perturbations from cosmolog- 
ical inflation, so we assume the Universe to be irrota- 
tional. As we are concerned about the present Universe, 
radiation is thus neglected. Moreover, the cosmological 
constant is also supposed to vanish, as we ask whether 
averaging could mimic a component of dark energy. Fol- 
lowing Buchert, we use physically comoving boundaries 
to thoroughly fix the averaging procedure. 

In the synchronous coordinates, the metric of the in- 
homogeneous and anisotropic Universe is ds 2 = — dt 2 + 
gij(t, x.)dx l dx : > . The spatial average of an observable 
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0(t, x) at time t is defined as 
1 



(0)d = y^rx J 0(t, x) x /det^dx. 



(1) 



Vo(t) = J D y/detgijdx is the volume of a comoving do- 
main D, introducing an effective scale factor 
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The subscript denotes the present time. The effective 
Hubble rate is thus defined as Hp = cld/cle, = (9)d/3 (9 
being the volume expansion rate) (T^ |. 

Effective Friedmann equations for a dust Universe fol- 
low from averaging Einstein's equations [l2l ] , 
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Here /j e ff and p e ff are the energy density and pressure of 
an effective fluid, 
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where p is the energy density of dust. (Q) p = §((# 2 ) d — 
Wl>)~ 2(c 2 ) p denotes the kinematical backreaction (a 2 
being the shear scalar) and (TV) p the averaged spatial 
curvature. They are related by an integrability condi- 
tion 
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We further define an effective equation of state, 

(TV) D -3(Q) D 
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So we find that cosmological backreaction gives rise to a 
nontrivial equation of state, even for a dust Universe [l6[ . 

Alternatively, we may map this effective fluid on a 
model with dust and dark energy. Let n be the num- 
ber density of dust particles, and m be their mass. For 
any comoving domain, (n)p = (n) p (aD /ap) 3 • In 
the dust Universe, p{t,x) = mn(t, x), and we identify 
Pm = (p)p = m(n)p. From Eq. (TJ|, dark energy is 
consequently pd c = —((Q)p + (7£)_d)/(167tG), with the 
relevant equation of state reading 
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It is — 1, iff (Q)p = — ^(R)p [1J|, corresponding to a 
cosmological constant A = (Q)p. 

Equations and ^ are not closed, as the four un- 
known variables (Q) p , (TV) p, (p)p and a jy are con- 
strained by only three equations. Below, we close these 



dynamical equations for the averaged Universe by means 
of cosmological perturbation theory. 

We wish to estimate the scale dependence of (Q)p, 
(TV)p, (p)pi Hp and w e g. We start from a spatially 
flat dust model. In the comoving synchronous gauge, 
the linear perturbed metric is ds 2 = — dt 2 + a 2 (t)[(l — 
2"$>)6ij + (didj — ^Sij A)x]dccMx J . Here, its scale factor 
a(t) (ao = 1) is different from old, and their relation 
was provided in Ref. ^ and % are the scalar metric 
perturbations, and A is the three-dimensional Laplace 
operator. The solutions for 'J and \ are given in terms 
of the time-independent peculiar gravitational potential 
p(x): * = ±A<^ /3 i 2 / 3 + | <p and x = -3wtt /3 t 2 / 3 (only 
growing modes are taken into account) j 1 61 ] - Moreover, 
<p is related to the hypersurface-invariant variable £ (l7j 

by C = iA^/V/i- |V. 

Following Ref. [16| , we use the metric perturbations at- 
tained from linear perturbation theory together with the 
non-perturbative integrability condition to obtain the av- 
eraged physical observables up to second order. We focus 
on the dominant contributions from the growing modes 
and neglect the decaying ones, since we are interested in 
the late time effects of cosmic averaging. Thus, we find 
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with B(tp) = (d l (d % (pAip) - d^djipd^dup)) - i(Aip} 2 and 
(O) = J D Odx/ J D dx. We see from Eqs. (ff|) - (|TT]l that 
these quantities are polynomials of surface terms. Thus, 
all information is encoded on the boundaries of the co- 
moving domain D. The temporal dependence of these 
averaged quantities can be found in Ref. and their 
leading terms are gauge-invariant (l6j . 

Our perturbative results suggest to write (Q) d and 
(TV) p in a Laurent series of an- (Recently, a power- 
law ansatz for the integrability condition was investi- 
gated in Ref. [l|.) We know from Eqs. © and © 
that (Q)p and (TV)p start from different powers: a^ 1 
and a D 2 , so (Q)p = E„=-i Qn(f£) n and (TV) D = 
5TJ n __ 2 The integrability condition then con- 

nects the coefficients: (n + 6)Q n + (n + 2)TZ n = 0. 
Thus, Qo — — at third order in perturbation the- 
ory. Therefore, cosmological backreaction can mimic a 
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cosmological constant, but induces extra terms as well. 
The third order results will be presented elsewhere . 

The effect of cosmological backreaction in the early 
Universe is tiny and is undistinguishable from that of a 
homogeneous curvature, as Wdc — ► — 1/3 when ao —* 0. 
This result seems inconsistent with our intuition of a van- 
ishing cosmological backreaction at early times, suggest- 
ing that should also vanish. However, as we have seen 
above, cosmological averaging gives rise to extra degrees 
of freedom in the dynamics of the averaged Universe. 

The effect of averaging over a typical domain is pro- 
vided by the ensemble average. From Eq. we find 
= x(^-) 2 *o( A( ^) 2 > °- This means that cosmo- 
logical backreaction is expected to lead to a positive def- 
inite equation of state. However, we should pay atten- 
tion that u^ff is of second ord er, but the root of its vari- 
ance [Varies)] 1 / 2 = f^i 2 l ("(A^)2) 1 /2 ig of first Qrder 

and therefore larger than w^s. Thus, the possibilities 
of w e g < and the effective acceleration of the aver- 
aged Universe cannot be easily excluded. Looking only 
at mean values of the ensemble obviously causes an un- 
derestimation of the possible backreaction effects, as we 
often observe just one particular domain in the Universe. 

We now turn to estimate the order of magnitude of 
cosmological backreaction as a function of the averaging 
scale : 



and is expected to fluctuate on small scales. Thus, 
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i.e., each derivative in B(ip) contributes a factor 1/r. 
Also, (ip 2 ) in the Fourier space is estimated as the power 
spectrum V v . Since (p is constant in time, and £ w — 5(^/3 
on superhorizon scales, we can identify today's V v with 
97^/25. Similar estimation works for the other two terms 
in B(ip). 

The kinematical backreaction induces 10% and larger 
modifications if \{Q)j}/4ttG{p)d\ ^ 0.1. This happens if 
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< - Mpc. 



(14) 



For observations at z -C 1, tq < 30 Mpc (h = 0.7). 

The averaged spatial curvature (IV) d is the most im- 
portant correction to energy density. The criterion for 
the scale, at which its effect emerges, is estimated analo- 
gously by 
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We show that cosmological averaging We find effects lar S er than 10% within 



produces important modifications to local physical ob- 
servables and determine the averaging scale, at which 
corrections show up at a 10% level. 

Effective acceleration of the averaged Universe occurs 
if pefi + 3p c ff < 0, i.e., (Q) D > 4irG(p) D . From Eqs. © 
and (©, we have 
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with i? H = 2.998 x lO 3 ^" 1 Mpc being the present Hub- 
ble distance. In Eq. (fl^j) . we can safely use the results 
for the background Unvcrsc: Od/cid = 1/(1 + z ) an d 
<o = 2_Rh/3, because B(ip) is of second order. Since the 
ratio in Eq. (I12p is dimensionlcss, a dimensional analysis 
immediately implies \ (Q)d/AtvG(p)d\ ^ (-Rh/?") 4 , where 
for an almost scale-invariant power spectrum the unique 
relevant scale is the averaging scale r. The order of mag- 
nitude of Eq. (fT2)) can be estimated as 



At small redshifts, rn < 77 Mpc. Furthermore, effects 
above 1% are expected up to a scale of ~ 240 Mpc. Note 
that the curvature of the Universe has been measured 
at the few per cent accuracy in the cosmic microwave 
background (CMB) [203. :t was shown in Ref. [H| that 
even small curvature might affect the analysis of high-z 
SNe significantly. 

Finally, we turn to the Hubble rate. To go beyond the 
order of magnitude estimates above, we calculate the en- 
semble mean and its variance (cosmic variance) of the rel- 
ative fluctuation of the Hubble rate 8h = (Ho — Hq)/Hq. 
Before doing so, let us stress that the analogous order 
of magnitude estimate for 5r agrees with the result for 
[Var((5_f/)] 1 / 2 given below up to a factor of ~ 2. For a 
spherical domain of radius r, we find from Eq. (TIT))) . 
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Vc, = 2.457 x 10~ 9 is the dimensionless power spec- 
trum [201 . We pick the second term in B(ip), 
(d l (djipd : 'diip)}, to demonstrate how to obtain this esti- 
mate. In the Fourier space, d l <p — ► ik l tp ~ (p/r. The lat- 
ter step comes from the observation that only structure 
of the size of the averaged volume cannot be averaged 
out. At much smaller scales, structures contribute a neg- 
ligible amount to B(tp), because it is not positive definite 
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with V — 47rr 3 /3 (a top-hat window). Above, we in- 
troduce the dimensionless power spectrum as </?ki Pk 2 = 
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2ir 2 6{k 1 +k 2 )7yfci)/fc? (fc= |k|). So 

Vai " (fe) = 144^ (1 + ^ {^j [ ^V v { X /r)Jl /2 { X ). 

(18) 

J 3 / 2 (a;) is the Bessel function of first kind (x = kr). For a 
scale-invariant power spectrum, we must introduce an ul- 
traviolet cutoff V v {k) — V ip e~ k / kc . No cutoff is required 
for a red-tilted spectrum V v {k) = V v (k/ko) nB ~ 1 (n s < 1 
being the spectrum index), consistent with WMAP5 [20| . 
Here, let us stress that although [Var (Sh)] 1 / 2 is only a 
hrst order quantity, the next contribution is already of 
third order, if we consult the perturbed metric to second 
order. Since we constrain our attention to the leading 
order effects, these higher order terms are negligible 

Now we can link the effect of cosmological backreaction 
in Buchert's setup (evaluated in a perturbative approach 
up to second order) to actual cosmological observations. 
The trick is to consider the scale dependence but not the 
time dependence. The value of the relative fluctuation of 
the Hubble rate in Eq. (jTTJ) is dominated by its variance, 
and thus the sign of the observed value of 5h cannot be 
predicted. A comparison of the mean and the root of the 
variance of Sh tells us that perturbation theory breaks 
down below ~ 20 Mpc. 

The scale dependence of the cosmic variance of Sh has 
previously been studied in the context of Newtonian cos- 
mology [22|, HH, largely based on CDM simulations. In 
this setting, the variance of Sh is due to peculiar motions 
(besides sampling variance and observational errors). In 
a relativistic and comoving approach, peculiar velocities 
vanish identically, and the cosmic variance of the Hubble 
rate turns into a curvature effect, because Eqs. ([5]) and 

G7> give Var(S H ) oc (K) 2 D . 

In Fig. (fTJ) , we compare the relativistic (correct up to 
second order) result Eq. (fT5)) to Newtonian "standard 
CDM" case in Ref. We find that up to - 400 Mpc, 
our results for scale-invariant power spectra (fc c = 1 /kpc 
corresponding to a typical cutoff in CDM simulations and 
1/pc to the physical cutoff in the primordial CDM spec- 
trum) agree with Newtonian simulations. This agree- 
ment is not unexpected, as metric perturbations and pe- 
culiar velocities are small at ~ 100 Mpc scales. 

The consistency between the relativistic and Newto- 
nian approaches encourages the comparison of our per- 
turbative results with experimental data. We compare 
Eq. (fTB"|) with observations from the HST Key Project @. 
We use 54 individual measurements of H in the CMB 
rest frame (corrected for local flow) from SN la and the 
Tully-Fisher relation (Tabs. (6) and (7) in Ref. H). We 
have checked explicitly that the SN and Tully-Fisher 
measurements of Hq are consistent with each other, 
while we cannot confirm that for the fundamental plane 
method and thus dropped them from a former analysis. 

We restrict our analysis to objects between 31.3 to 
467.0 Mpc, as Eq. flU) can be trusted only above 30 Mpc. 
Be r,, Hi and er^ the distance, Hubble rate and la er- 
ror for the i'th datum, with distances increasing. We 
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FIG. 1: Scale dependence of the cosmic variance of the Hubble 
rate. Data are from the Newtonian CDM model in Ref. [23I ] . 
with h = 0.5, £lm = 1 and a COBE-normalized power spec- 
trum. Thick and thin lines correspond to the relativistic result 
Eq. (|18[1 for a scale-invariant power spectrum with cutoffs at 
fc c =l/kpc (simulation) and 1/pc (physical), respectively. 
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FIG. 2: Relative fluctuation of the Hubble rate from cosmo- 
logical backreaction and its cosmic variance band (thick lines) 
compared to the empirical mean and variance of 5h obtained 
from the HST Key Project data |5j as a function of averaging 
radius. The thin line shows the ensemble mean of Sh- The 
band enclosed by the thick lines indicates the effect of the 
inhomogeneities (oc 1/r 2 ), and the dashed lines are the effect 
from sampling with given measurement errors in an otherwise 
perfectly homogeneous Universe. 



calculate the mean distance for the nearest k objects 

by f k = £) i= i X),=i ft) witn weights g t = H$/af. 
An analogue holds for the averaged Hubble rate Hk, i.e., 
Hd for different subsets. The empirical variance of each 
subset is -a\ = - H k ) 2 ]/[Hi(k - 1) £* =1 9i ]. 

Notice that Eq. (fT5|) is insensitive to global calibration 
issues. 

The comparison of the result Eq. (fT5|) with the HST 
Key Project data is shown in Fig. We now normalize 
to the WMAP5 best-fit power-law spectrum, with pivot 
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fco=0.002/Mpc and spectral index n s = 0.960 and use 
H = 72 km/s/Mpc [H. We see that the theoretical 
band matches the experimental data well, without any fit 
parameter in the panel. Moreover, we see from Fig. @ 
that the value of Sh is positive within ~ 100 Mpc. This 
is consistent with the result in a recent paper [24|] that we 
are located in a 200 - 300 Mpc underdense void, which 
is expanding faster than the global Hubble rate. 

Before we can claim that we have observed the ex- 
pected 1/r 2 behavior in Eq. (fT8"|) and thus the evidence 
for cosmological backreaction, we must make sure that 
statistical noise cannot account for it. In the case of a 
perfectly homogeneous coverage of the averaged domain 
with standard candles, we would expect a 1/r 3 / 2 behav- 
ior. In Fig. ([2|), we show the statistical noise for the ac- 
tual data set (1/Q3j=i sO 1 ^ 2 )' which is smaller than our 
result Eq. (fT%|) . It turns out that the sampling noise for 
this small data set is still too large to claim that the inho- 
mogeneity of the Universe can be detected in the relative 
fluctuation of the Hubble rate observed by the HST Key 
Project. However, it is fully consistent with our theoret- 
ical expectations. Actually the fluctuation 5h appears 
to be smaller than expected, and one might wonder why 
that is so. From the theoretical expectation plotted in 
Fig. (|2l), we find that at ~ 40 (60) Mpc, the value of H D 
differs from its global value 72 km/s/Mpc (WMAP5) by 
about 10% (5%), whereas the expected variance for a per- 
fectly homogeneous and isotropic Universe is 8% (2%). 

A similar analysis of the Hubble diagram was pioneered 
in Refs. (2f| [26|, [27], [28j], in which the velocity field of the 
local Universe and its influence on the correlated fluctua- 
tions in luminosity distance and the Hubble rate was ana- 
lyzed. Two essential differences to this work are that our 
analysis includes effects to higher orders and we study 



the scale dependence of the averaged observables. Al- 
though the relative fluctuation of the Hubble rate was 
not explicitly analyzed in Refs. 0, [H, it seems to 
us that our results are consistent with those findings. 

To summarize, we argue that cosmological averaging 
(backreaction) gives rise to observable effects up to scales 
of ~ 200 Mpc. However, it is not sufficient to explain the 
observed accelerated expansion at this point. 

We find a hierarchy of backreaction effects. The av- 
eraged spatial curvature (11) d leads to 10% (1%) effects 
up to ~ 80 (240) Mpc in a dust model with h = 0.7. 
Below ~ 40 Mpc, the cosmic variance of the Hubble 
rate is larger than 10%, which coincides with the esti- 
mate from the effect of peculiar motions in Newtonian 
setup. Within ~ 30 Mpc, the kinematical backreaction 
(Q)d-, due to second order perturbations caused by lo- 
cal inhomogeneities and anisotropies, enters the game. 
Cosmological backreaction may put some of the steps on 
the cosmological distance ladder in question, as they are 
deeply in the domain of large backreaction, i.e., a large 
fluctuation between small averaged volumes. 

Our findings call for revisiting local observations, like 
galaxy redshift surveys, in terms of possible backreaction 
signatures. The large scale physics of primordial CMB 
anisotropies is not affected. However, this statement can- 
not be made for secondary effects, e.g., the late integrated 
Sachs- Wolfe effect. 
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